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Abstract. We investigate ways of accurately simulating the propagation of ener- 
getic charged particles over small times where the standard Monte Carlo approx- 
imation to diffusive transport breaks down. We find that a small-angle scattering 
procedure with appropriately chosen step-lengths and scattering angles gives ac- 
curate results, and we apply this to the simulation of propagation upstream in 
relativistic shock acceleration. 
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1. Introduction 

Relativistic charged particle transport in magnetized astro-physical 
plasma is strongly affected by magnetic irregularities, and may be 
approximated by diffusion. Diffusive transport of particles having speed 
v can be simulated by a three-dimensional random walk with steps sam- 
pled from an exponential distribution with mean free path A = 3D/v, 
where D (cm 2 s _1 ) is the spatial diffusion coefficient, followed by large- 
angle (isotropic) scattering after each step (e.g. Chandrasekhar 1943), 
and this gives good results for distances much larger than A. 

In diffusive shock acceleration at relativistic shocks problems arise 
when simulating particle motion upstream of the shock because the 
particle speeds, v, and the shock speed v s hock = c(l — l/7shock) 1//2 are 
both close to c, and so very small deflections are sufficient to cause 
a particle to re-cross the shock. Clearly, Monte Carlo simulation by 
a random walk with mean free path A and large-angle scattering is 
inappropriate here, and in Monte Carlo simulations of relativistic shock 
acceleration at parallel shocks Achterberg et al. (2001) consider instead 
the diffusion of a particle's direction for a given angular diffusion coeffi- 
cient Dq (rad 2 s _1 ). Similarly, for a given spatial diffusion coefficient D, 
Protheroe (2001) and Meli & Quenby (2001) adopt a random walk with 
a smaller mean free path, f C A, followed by scattering at each step 
by a small angle with mean deflection , 9 < l/7 s hock- See Bednarz 
Ostrowski (2001) for a recent review of relativistic shock acceleration. 
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2. Small-angle scattering and diffusion 

We consider propagation by small steps sampled from an exponential 
distribution with mean t <C A, followed at each step by scattering 
by a small angle sampled from an exponential distribution with mean 
9 <C 7r. The change in direction {61,62) may then be described as two- 
dimensional diffusion with angular diffusion coefficient Dg = Ovg/2 
(rad 2 s" 1 ) where vg = 6/i, and t = i/v such that Dg = 9 2 v/{2i). The 
time tiso, which gives rise to a deflection equivalent to a large-angle 
(isotropic) scattering, is determined by (<r| + c 2 2 ) 1//2 = \/4:Dgti so ~ 
7r/2, giving A ~ vt- lso oc £/# 2 and a spatial diffusion coefficient D oc 
£t>/# 2 oc v 2 /Dg. By using a Monte Carlo method it is straightforward 
to test this, determine the constant of proportionality, and thereby 
make the connection between diffusion and small angle scattering. 

The solution of the diffusion equation for a delta-function source 
in position and time q(f,t) = 5{r)5{t) and an infinite diffusive medium 
is a three-dimensional gaussian with standard deviation a = \/2Dt 
(Chandrasekhar 1943). The results from several Monte Carlo random 
walk simulations are shown in Fig. 1, from which we find that the 
expected dependence occurs for 9 < 5° at times t > 10 5 £/v. For this 
case, we see from Fig. 1 that a 2 — ► 2tv£/39 2 , and so we obtain the 
connection between small-angle scattering and diffusion theory, namely 
D pa £v/{39 2 ) pa v 2 /(QDg). 

3. Application to relativistic shock acceleration 

As viewed in the frame of reference of the upstream plasma, ultra- 
relativistic particles are only able to cross the shock from downstream 
to upstream if the angle 9 between their direction and the shock nor- 
mal pointing upstream is 9 < sin _1 (l/7 s hock), where 7 s hock = (1 — 
^s 2 hock)~ 1/2 and Ashock = fshock/c For highly relativistic shocks these 
particles cross the shock from downstream to upstream travelling al- 
most parallel to the shock normal. Similarly, having crossed the shock, 
only a very slight angular deflection, by ~ l/7 s hock is sufficient to return 
them downstream of the shock. This change in particle direction gives 
rise to a change in particle energy E' and momentum p', measured in 
the downstream plasma frame (primed coordinates), of 

K+l „ Pn+l _ 1-Pl2 COS gn+l ^ 

E' n ~ p' n l-/? 12 cos0 n 1 ' 

in an acceleration cycle (downstream — > upstream — > downstream), 
where j3\2 is the speed of the upstream plasma as viewed from the 
downstream frame. 
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Figure 1. a vs. time for a 3D random walk with isotropic injection at the origin 
at t — 0. Step-lengths £ were sampled from an exponential distribution with mean 
I followed by small-angle scattering with scattering angle 8 sampled from an expo- 
nential distribution with mean 9 (the numbers attached to the curves). Dashed line 
is a 2 = 2tv£/38 2 . Curves for 5°-20° result from 10 4 simulations; 4° curve results 
from 8 x 10 4 simulations (width shows statistical error). 

In "parallel shocks" the magnetic field is parallel to the shock 
normal, and so the the pitch angle Y> is the angle to the shock normal 
and v cos ip gives the component of velocity parallel to the shock. Thus 
the small-angle scattering method described in the previous section is 
used here to simulate particle motion upstream of a parallel relativistic 
shock, including the effects of pitch-angle scattering, for a given diffu- 
sion coefficient. We inject ultra-relativistic particles at the shock with 
downstream-frame energy E' travelling upstream parallel to the shock 
normal, i.e. 9q = 0. We follow a particle's trajectory until the shock 
catches up with it, and it crosses from upstream to downstream with 
an upstream-frame angle 9\ to the shock normal and a downstream- 
frame energy E[. The simulation was performed for 7 sn0 ck = 10, and 
five different mean scattering angles 9, to determine the maximum 9- 
value that can safely be used for accurate simulation. The resulting 
distributions of cos# and log(E[/E' ) are shown in Fig. 2, and show 
that in this application one requires 9 < 0.1/7 snO ck- Our results are 
quite consistent with those of Achterberg et al. (2001), who used a 
diffusive angular step A# s t < 0.1/7 s h oc k- 
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a. cos Oi distribution b. \og(E[/E' ) distribution 

Figure 2. Small-angle scattering simulation of excursion upstream in diffusive shock 
acceleration at a parallel relativistic shock with 7 s hock = 10. Results are shown 
for 10 5 injected particles and 9 — 10~ 2 /7 s h oc k (top histogram), 3 x 10~ 2 /7 s h O ck, 
10 _1 /7shock, 0.3/7 s hock, l/7shock and 3/7 s hock (bottom histogram). Note that the 
top three histograms are almost indistinguishable. 



4. Conclusion 

The standard Monte Carlo random walk approach to simulation of 
energetic charged particle propagation for a given spatial diffusion co- 
efficient D can be extended to apply accurately to times much less than 
\/v = 3D/v 2 by using a small-angle scattering procedure with steps 
sampled from an exponential distribution with mean free path £ = 9 2 X 
followed at each step by scattering with angular steps sampled from an 
exponential distribution with mean scattering angle 9 < 0.09 rad (5°). 
The spatial and angular diffusion coefficients are then D rs £v/(39 2 ) 
and Dq = 9 2 v/(2£), and are related by D « v 2 /(6De). In simulation of 
upstream propagation in relativistic shock acceleration one must use 
9 < 0.1/7 s hock to obtain accurate results. 
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